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ABSTRACT 

A  manual  flight  control  system  for  the  lateral 
directional  dynamics  of  a  modern  fighter  aircraft 
incorporating  thrust  vectoring  is  presented.  Design 
goals  are  posed  in  terms  of  maintaining  acceptable 
flying  qualities  during  high  a  maneuvering  while  also 
achieving  robustness  to  model  parameter  variations  and 
unmodeled  dynamics  over  the  entire  flight  envelope. 
The  need  for  gain  scheduling  is  eliminated  by  using  an 
inner  loop  dynamic  inversion/outer  loop  structured 
singular  value  (|Li)  -  synthesis  control  structure  which 
separately  addresses  operating  envelope  variations  and 
robustness  concerns,  respectively.  Performance 
objectives  are  based  on  commanding  sideslip  angle  and 
stability  axis  roll  rate.  Realistic  representations  of  both 
structured  (real  parametric)  and  unstructured  uncertainty 
are  included  in  the  design/analysis  process.  A  flight 
condition  dependent  control  selector  maps  generalized 
controls  to  physical  control  deflections,  considering 
actuator  redundancy,  effectiveness,  and  saturation  issues. 
An  angle  of  attack  dependent  command  prefilter  shapes 
commands  to  produce  desired  responses.  Structured 
singular  value  analysis,  low-order  equivalent  system 
(LOES)  fits,  and  linear  step  responses  demonstrate 
satisfaction  of  design  goals.  Simulation  shows 
excellent  control  at  both  low  and  high  angles  of  attack 
without  gain  scheduling. 

INTRODUCTION 

Although  many  of  today's  fighters  can  achieve 
a  'first  shot’  on  the  sheer  ability  to  acquire  a  target  from 
many  miles,  failures  and  rules  of  engagement  make 
superior  maneuverability  a  primary  mode  to  gain 
advantage.  With  the  advent  of  unconventional  control 
effectors  like  thrust  vectoring  and  forebody  vortex  flow 
control,  modern  aircraft  will  be  able  to  gain  this 
advantage  by  maneuvering  at  and  beyond  the  stall  angle 
of  attack,  giving  pilots  superior  survivability,  and 
increased  first-shot  opportunities^.  This  potential  for 
high  angle  of  attack  flight  presents  many  challenges  to 
the  control  designer.  The  control  law  must  incorporate 
both  the  conventional  and  unconventional  effectors, 
provide  precise  control  and  handling  qualities,  and 
provide  robust  stability  in  the  face  of  nonlinearities 
(saturation,  inertial  coupling)  and  increased  uncertainty 
(lack  of  data ). 


*  Stability  and  Control  Engineer,  Member  AIAA 

**  Captain,  U.S.  Air  Force 
t  Aerospace  Engineer,  Associate  Fellow  AIAA 
This  is  declared  a  work  of  the  U.S.  Government  and  is 
not  subject  to  copyright  protection  in  the  United  States. 


In  recent  years,  modern  multivariable  methods 
have  been  applied  to  a  variety  of  model  aircraft  which 
had  the  physical  ability  to  fly  at  high  angles  of  attack. 
A  robust  Hoo  controller  in  an  inner/outer  loop  structure 
has  been  designed  for  a  thrust  vectored  aircraft  at  a 
single  flight  condition^’^.  A  Herbst-like  maneuver  was 
done  to  demonstrate  the  maneuverability.  Sparks'^ 
applied  |i-synthesis  to  the  same  aircraft,  incorporating 
flying  qualities  and  accounting  for  structured 
uncertainties.  Along  those  lines,  a  Wright  Laboratory 
technical  report^  documents  the  design  of  a  |i-synthesis 
controller  for  full  envelope  robust  stability  and 
performance  for  high  angle  of  attack  flight  using  the 
inner/outer  loop  control  structure.  Buffington  et  al^ 
designed  dynamic  inversion/|a-synthesis  inner/outer  loop 
control  for  a  thrust  vectored  F-18.  They  also 
incorporated  a  control  allocation  scheme  to  minimize 
adverse  control  power  saturation  effects  and 
demonstrated  performance  through  a  Herbst-like 
maneuver.  NASA  Langley^’^  has  done  work  on  feed 
forward  control  architecture  for  the  F-18  HARV  and  did 
piloted  simulation  to  evaluate  high  angle  of  attack 
handling  qualities. 

The  use  of  thrust  vectoring  has  not  been  the 
only  form  of  performance  enhancement.  Adams  et  alP 
developed  a  nonlinear  control  strategy  for  the  control  of 
forebody  vortex  flow  for  the  F-16  through  the  use  of  a 
forebody  blowing  system.  Results  showed  that  high 
angle  of  attack  performance  was  significantly  improved. 

As  previously  mentioned,  the  high  angle  of 
attack  region  is  very  nonlinear  compared  to  the  low  to 
mid-level  angle  of  attack  regions  and  has  precipitated 
research  in  the  area  of  nonlinear  controls.  Specifically, 
nonlinear  dynamic  inversion  has  been  applied  to  the 
F-18  HARV^^’^^  for  control  over  a  wide,  high  angle  of 
attack  envelope.  A  dynamic  inversion  approach  has 
also  been  used  to  develop  control  laws  to  control  the 
X-29  at  high  angles  of  attack^ 

Control  law  design  has  not  been  the  only  area 
of  nonlinear  research.  Another  major  issue  of  high 
angle  of  attack  flight  is  control  saturation.  Various 
methods  to  prevent  instability  due  to  saturation  have 
been  examined.  Bugajski  and  Enns^^  examined  a 
method  in  which  loop  bandwidths  are  reduced  based  on 
an  achievable  subspace.  Similarly,  Durham proposed 
a  method  of  allocating  control  effectors  such  that  the 
maximum  moment  is  achieved  given  an  achievable 
moment  subspace.  Lateral  command  scaling,  and 
longitudinal  prioritization  were  examined  by  Buffington 
et  al^.  All  of  these  methods  were  shown  to  have 
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satisfactory  results.  Also,  research  done  by  Hammett  et 
alM  has  found  that  the  method  of  daisy-chaining  can 
produce  excessive  phase  lag  in  the  longitudinal  control 
of  the  VISTA/MATV  F-16,  reducing  its  stability  and 
performance  envelope. 

The  main  contribution  of  this  paper  is  the 
integration  of  certain  aspects  of  the  previously 
mentioned  research  into  application  to  a  main-line 
fighter  platform  that  has  been  modified  and  flight  tested. 
This  paper  details  application  of  a  modified  version  of 
the  three  part  design  approach  of  Adams  et  al^  to 
manual  flight  control  system  design  for  a  model  of  the 
lateral-directional  dynamics  of  the  VISTA/MATV  F-16 
flight  test  aircraft. 


Controller 


Figure  1  Modified  Three  Part  Control  Structure 

The  structure  (Figure  1)  consists  of  an  inner  loop  design 
that  equalizes  the  dynamics  of  the  plant  over  the 
envelope  of  concern,  an  outer  loop  controller  that 
conveniently  addresses  performance  and  robustness 
issues,  a  command  prefilter  for  handling  qualities,  and  a 
control  selector  that  implements  generalized  controls, 
allowing  the  designer  to  prioritize  and  combine  control 
surfaces  as  desired.  The  inner  loop  control  is  a  modified 
form  of  dynamic  inversion  consisting  of  a  gain 
matrix  that  is  flight  condition  dependent,  and  the  outer 
loop  controller  is  a  single  |a-synthesis  controller  found 
at  a  single  design  point.  A  major  benefit  of  this 
structure  is  that  the  inner  loop  equalization  is 
accomplished  through  a  simple  table  lookup  of 
aerodynamic  data  and  allows  a  single  dynamic  outer 
controller  to  work  over  the  entire  flight  envelope, 
avoiding  the  difficulties  associated  with  gain  scheduling. 
In  addition,  the  |i-synthesis  outer  loop  can  address  the 
performance  and  robustness  issues  associated  with  the 
problem.  The  control  selector  employs  the  concept  of 
pseudo-controls  to  facilitate  efficient  use  and 
combination  of  control  power,  and  is  found  through 
matrix  operations  involving  the  control  effectiveness 
matrix  and  designer  chosen  weights.  The  control 
selector  also  structures  the  inner  loop  so  that  dynamic 
inversion  is  readily  implemented.  The  prefilter  is  used 
to  provide  the  desired  handling  qualities  characteristics 
by  shaping  the  stick  commands  based  on  the  dynamics 
provided  by  the  outer  loop  controller.  Since  the  desired 
response  of  the  aircraft  is  different  for  different  angles  of 
attack,  the  parameters  of  the  prefilter  are  varied  with 
angle  of  attack. 


In  the  following  sections,  a  description  of  the 
VISTA  aircraft  model  is  given  followed  by  design 
objectives.  The  control  architecture  and  design  are 
described  followed  by  linear  robustness,  handling 
qualities,  and  time  domain  analysis. 

AIRCRAFT  MODEL 

Wright  Laboratory's  Variable  Stability  In- 
Flight  Simulator  Test  Aircraft  (VISTA)  is  a  modified 
F-16  with  the  capability  to  simulate  and  test  advanced 
aircraft  configurations  and  flight  control  concepts. 
Recently,  the  VISTA  F-16  was  modified  by  the  addition 
of  a  multi-axis  thrust  vectoring  (MATV)  nozzle  (Figure 
2)  to  facilitate  investigation  of  the  tactical  utility  of 
high  a  flight.  The  primary  aerodynamic  control 
surfaces  are  left  and  right  horizontal  stabilators,  left  and 
right  flaperons,  and  the  rudder.  The  thrust  vectoring  can 
be  controlled  in  pitch  and  in  yaw  to  provide  additional 
control  power. 


Figure  2  VISTA/MATV  Aircraft 


A  nonlinear  simulation  model  exists  as 
modular  FORTRAN  code.  The  various  modules 
describe  the  atmosphere,  nonlinear  equations  of  motion, 
aerodynamics,  engines,  sensors,  thrust  vectoring, 
control  system,  and  actuators.  The  simulation  can 
generate  a  ten  state  linear  model  at  any  flight  condition 
with  any  given  set  of  inputs  and  outputs.  This  ten  state 
model  is  then  broken  down  into  the  following  three 
state  lateral-directional  state  space 


Y p  sin  a  -  cos  a 


^DE 
L  ^DE 


^DF 

L  (j  L  jr 

^DF 


«DF 


"R 


^YTV 

^YTV 


^DE 

^DF  1 

/ytv. 

(2.1) 


where  (3  is  yaw  angle,  p  is  body  axis  roll  rate,  r  is  body 
axis  yaw  rate,  5de  is  differential  elevator  deflection, 
8dF  is  differential  flaperon  deflection,  5r  is  rudder 
deflection,  and  6ytV  is  yaw  thrust  vectoring  deflection. 
Side  force,  roll  moment,  and  yaw  moment  are 
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represented  by  Y,  L,  and  N  respectively.  All  three 
variables  are  assumed  to  be  available  for  feedback  and  all 
angles  are  in  degrees  or  degrees/second.  One  weight  and 
store  configuration,  corresponding  to  two  AIM-9 
missiles  on  the  wing  tips,  and  trim  throttle  setting  are 
assumed. 

PERFORMANCE  OBJECTIVES 

The  performance  objectives  for  this  problem 
involve  the  robust  tracking  of  sideslip  and  velocity 
vector  roll  rate  command  with  desired  handling  qualities 
specified  in  MIL-STD-1797A.  Although  the 
requirements  given  in  the  standard  are  not  valid  for  high 
angle  of  attack  flight,  enough  guidance  was  given  for 
Buffington  et  al.^  to  extrapolate  to  high  angle  of  attack 
estimates.  Evaluation  of  flying  qualities  is  performed 
by  obtaining  a  Low  Order  Equivalent  System  (LOES) 
fit  of  the  form  in  (3.1)  to  the  closed  loop  system. 

^lat  (s+^)(s+^)(s^  +2Cj(b^s+(Uj2  ) 

^  r  ^  s 

In  this  equation,  /i  is  stability  axis  roll  rate,  5iat  is 
command  input,  Ti-  is  roll  mode  time  constant,  Tg  is 
spiral  mode  time  constant,  is  dutch  roll  damping,  coj 
is  dutch  roll  frequency,  and  t[x  is  the  equivalent  time 
delay.  The  other  values  in  (3.1)  are  gains  or  zero 
dynamics  which  are  not  specified  in  the  standard. 
Calculating  the  values  for  (3.1)  involves  a  three  step 
process  that  is  presented  in  Adams  et  al.^  and  is  as 
follows.  First,  the  transfer  function  from  lateral 
command  to  stability  axis  roll  rate  is  fit  to  the  form  of 
equation  (3.2)  to  find 


V 

p 

(3.4) 

fd®j>0,4 

(3.5) 

m ,  >1.0  rad  /  s 
d 

(3.6) 

spiral  time  to  double  <  12  sec (3.7) 

Tj.  <1.0  second 

(3.8) 

Time  delay  <0.1  second 

(3.9) 

Since  it  is  desired  to  have  different  handling 
qualities  at  various  angles  of  attack,  these  values  have 
been  relaxed  for  some  flight  conditions.  In  particular, 
the  roll  time  constant  is  relaxed  to  1.4  seconds  (level  2) 
at  high  angles  of  attack.  This  makes  physical  sense 
since  there  is  less  control  power  at  high  angle  of  attack 
making  the  responses  slower.  Also,  due  to  the 
susceptibility  to  depart  at  high  angles  of  attack,  the 
dutch  roll  damping  must  be  greater  than  0.8,  and  the 
dutch  roll  frequency  will  lie  near  the  edge  of  level  one 
requirements. 

CONTROL  SELECTOR  DESIGN 

The  control  selector  accomplishes  two 
important  things  in  the  three  part  control  structure.  The 
first  is  that  it  takes  generalized  rate  commands  and 
transforms  them  into  actuator  deflection  commands. 
The  second  is  that  it  allows  for  the  redistribution  of 
control  power  between  the  various  effectors  without 
affecting  the  linear  closed  loop  performance.  As  a 
result,  the  designer  can  then  combine  and  prioritize 
control  usage  based  on  effectiveness  or  limits. 

The  control  selector  is  found  by  redefining  the 
control  contribution  to  the  state  equations^, 

B5  =B*6*.  (4.1) 


/i(s)  ^ 

(s+^) 


(3.2) 


The  next  step  is  to  use  the  high  order  transfer  function 
of  rudder  pedal  to  sideslip  to  find  the  dutch  roll  values  in 

(3.3). 


B  and  5  are  the  actual  control  effectiveness  matrix  and 
control  vector,  and  B*  and  5*  are  the  generalized 
equivalents.  The  actual  control  can  be  defined  in  terms 
of  the  generalized  control  by  introducing  a 
transformation  matrix,  T,  such  that 

8  =  T5*.  (4.2) 


-TgS 

Pis)  ^ 


(3.3) 


We  call  T  the  control  selector.  It  may  be  calculated  as 
T  =  N(BN)^B*.  (4-3) 


The  roll  constant  and  dutch  roll  values  from  (3.2)  and 

(3.3)  are  substituted  into  (3.1)  and  fixed.  Another  fit  is 
then  done  to  find  the  equivalent  time  delay,  numerator 
roots,  and  spiral  mode  time  constant.  The  ranges  for 
the  various  roll  mode  values  that  produce  level  one 
flying  qualities  for  a  fighter  type  aircraft  are  shown  in 

(3.4)  to  (3.9) 


The  operation  ( )^  is  a  pseudo-inverse  and  N  is  a  matrix 
that  may  be  used  to  combine  controls  or  prioritize 
individual  control  channels  in  the  case  of  redundant 
effectors.  Because  the  B  matrix  in  (2.1)  is  a  function  of 
flight  condition,  the  control  selector  is  a  function  of 
Mach,  altitude,  and  angle  of  attack. 

The  generalized  inputs  and  actual  controls  for 
this  design  are: 
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(4.4) 


The  actual  control  effectiveness  matrix  is  given  in  2.1 
and  2.2,  and  the  generalized  control  effectiveness  matrix 
is  the  three  by  two  matrix: 


B 


* 


'0  0” 
:  1  0 
0  1 


(4.5) 


The  origin  of  this  control  effectiveness  matrix  is  a  direct 
result  of  the  dynamic  inversion  and  is  discussed  in  the 
next  section. 

Specification  of  the  control  weighting/ganging 
matrix  N  is  now  considered.  The  driving  force  behind 
the  choice  of  N  for  the  lateral-directional  case  is  the 
actuator  dynamics  and  limits  which  are  shown  in  (4.6), 
(4.7),  and  in  Table  I. 

Aerodynamic  Effectors: 


since  consideration  must  be  given  to  the  fact  that  the 
tails  are  more  important  for  keeping  longitudinal 
stability  and  performance.  In  the  case  of  the  yawing 
moment  generators,  it  can  be  seen  from  Table  I  that  the 
rate  and  position  limits  for  the  rudder  are  roughly  twice 
that  of  the  thrust  vectoring.  Since  it  is  desired  that 
these  two  effectors  limit  at  around  the  same  point,  the 
thrust  vectoring  command  should  be  half  of  the  rudder. 
With  these  things  in  mind,  the  N  matrix  is 


ro.25 


0 


0  " 
0 

1.0 

0.5 


(4.8) 


An  example  is  shown  for  angle  of  attack^  6  degrees, 
altitude=20000  ft,  and  Mach=0.5.  The  B  matrix  is 


"  0.0187 
-15.0534 
-1.6781 


0.0019 

-20.5721 

-0.3163 


0.0204 

4.4799 

-1.8053 


0.0061  " 
-0.0890 
^.6049 


(4.9) 


_ (20.2)(144.8)(71.4)^ _ 

(s+20.2)(s+144.8)(s^+2(0.736)(71.4)s+(71.4)^) 

Thrust  Vectoring  Nozzle: 

_ (21.3)0 18.8)(26.3)^ _ 

(s+21.3)(s+118.8)(s^+2(0.444)(26.3)s+(26.3)2) 


(4.6) 


(4.7) 


Table  I  Actuator  Limiting  Cl 

laracteristics 

Control 

Effector 

Rate  Limit 
(Deg/s) 

Position  Limit 
(Deg) 

Horizontal 

Tails 

±  50 

±  21 

Flaperons 

±  52 

+20,-23 

Rudder 

±  120 

+  30 

Pitch  and  Yaw 
Nozzle 

±  45 

+  17 

In  early  designs,  N  was  chosen  to  be  a  four  by  four 
diagonal  matrix  with  weightings  on  the  diagonal 
corresponding  to  the  different  actuators.  This  was  found 
to  give  the  system  too  much  freedom  to  use  the 
actuators.  As  a  result,  the  tails  and  flaperons  would 
work  against  each  other  but  still  produce  good  linear 
results.  As  soon  as  actuator  limits  were  introduced,  the 
system  would  go  unstable.  Since  the  primary  modes  of 
control  are  the  rolling  and  the  yawing  of  the  aircraft,  the 
actuators  are  ganged  into  two  sets.  The  tails  and 
flaperons  are  ganged  since  they  primarily  produce  body 
axis  rolling  moments,  and  the  rudder  and  thrust 
vectoring  are  ganged  since  they  both  produce  primarily 
yawing  moments.  The  tails  are  weighted  so  that  they 
contribute  one  quarter  the  control  power  of  the  flaperons 


which  when  used  in  (4.3)  produces 


"-0.0097  -0.0203" 
-0.0386  -0.0813 
0.0135  -0.4460 

0.0067  -0.2230 


(4.10) 


One  can  see  that  the  primary  roll  effector  is  the 
flaperons,  and  the  primary  yaw  effector  is  the  rudder. 
Also,  the  differential  tail  is  commanded  one  fourth  the 
amount  of  the  flaperons,  and  the  thrust  vectoring  is 
commanded  one  half  the  rudder  command  for  any  given 
input  of  p  and  f . 

INNER  LOOP  DESIGN 

The  purpose  of  the  inner  loop  is  to  equalize  the 
plant  dynamics  over  the  entire  envelope  so  that  a  single 
dynamic  controller  can  be  used  in  the  outer  loop.  A 
convenient  way  to  do  this  is  through  dynamic  inversion 
which  allows  the  designer  to  specify  the  closed  loop 
linear  dynamics.  Although  nonlinear  effects  such  as 
inertial  coupling  can  be  included  in  the  formulation, 
this  work  only  dealt  with  the  nonlinearities  associated 
with  the  aerodynamic  parameters  which  vary  with  Mach 
number,  altitude,  and  angle  of  attack.  Neglecting  these 
nonlinearities  does  introduce  error,  but  since  our  linear 
simulation  does  not  include  the  coupling  in  the  state 
equations,  the  effects  will  not  be  seen  here.  Should 
these  effects  prove  to  be  significant  in  full  nonlinear 
simulations,  they  can  be  canceled  by  simply  adding  the 
coupling  terms  to  the  inversion  equations. 

Traditional  application  of  dynamic  inversion 
involves  choosing  the  same  number  of  controlled 
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outputs  as  there  are  independent  controls.  The  two 
generalized  controls  have  been  chosen  to  be  p  and  r . 
Therefore,  only  two  outputs  can  be  controlled  and  they 
are  chosen  to  be  p  and  r.  These  variables  must  only  be 
differentiated  once  for  the  controls  to  appear  in  the 
output  equations. 


sin  a 


-cos  Of 


‘0  O' 
I  0 
0  1 


(5.1) 


Since  we  can  only  control  two  outputs,  the  row 
associated  with  sideslip  is  neglected  resulting  in  (5.2). 


{*}-K "Jl'k  ?]{<:} 


(5.2) 


The  inverse  dynamics  control  law  for  the  inner  loop  can 
be  written  as 


A- 


-0.1412 

-21.378 

4.4460 


0.1052  -0.9925" 

-1.9647  0.2127 

-0.0712  -0.3288 


(5.7) 


The  desired  poles  for  the  system  are  -3.5+3.5i,  -3.5- 
3.5i,  and  -5.0,  and  the  eigenstructure  desired  is  to  have 
the  complex  poles  associated  with  the  dutch  roll  and  the 
real  pole  with  the  roll  mode.  The  poles  are  chosen  so 
that  the  inner  loop  dynamics  are  faster  than  the  desired 
outer  loop  dynamics  but  slow  enough  to  prevent 
excessive  control  usage.  The  feedback  gain  matrix  is 
found  to  be 


pl.l642 

-3.0419 

-0.1432' 

tl9.3287 

0.5740 

-6.2466^ 

(5.8) 


which  produces  the  following  closed  loop  dynamics: 


-0.3527  0.0987 

-0.2138  -5.0066 
23.7723  0.5027 


-0.9237" 

0.0695 

-6.5746 


(5.9) 


ftH: 

V  y 

where  v  represents  the  desired  closed  loop  linear 
dynamics,  and  the  subscript  ci  denotes  inner  loop 
commands.  Our  desired  dynamics  can  be  represented  by 
a  matrix  multiplying  the  state  vector  as  follows: 


Lp  L, 
Np  Nr 


(5.3) 


The  poles  associated  with  (5.9)  are  -3.466-H3.5015i, 
-3.466-3.5015i,  -5.002.  The  matrix  A(j  comes  directly 
from  this  by  taking  the  last  two  rows  making  our 
desired  dynamics 

*  1^1  \-O.2m  -5.0066  0.0695  1  rl  /r- 

[23.7723  0.5027  -6.5746 ^  ^ 


v=A^ 


^dl2 

^d22 


(5.4) 


The  inner  equalization  loop  can  then  be  represented  as  a 
linear  state  feedback  compensator  of  the  form: 


(5.5) 


where 


K 


eq 


^d2l~^P 


^dl2 

1 

< 

^d22“^p 

^d23~^: 

(5.6) 


Neglecting  the  sideslip  dynamics  does  introduce  some 
error  into  the  formulation,  but,  more  importantly,  it 
raises  an  issue  of  stability.  Since  the  number  of 
outputs  controlled  is  less  than  the  number  of  states  in 
the  plant,  stability  is  not  guaranteed.  Care  must  be 
taken  to  check  the  internal  stability  of  the  dynamic 
inversion  at  various  operating  points  to  insure 
stability^ 

The  actual  implementation  of  the  dynamic 
inversion  is  a  simple  process.  The  aerodynamic 
derivatives  needed  for  (5.6)  are  stored  in  a  tabular 
database  and  a  table  lookup  and  subtraction  are  done  to 
find  Keq.  The  inner  loop  generalized  commands  are 
then  obtained  by  multiplying  the  feedback 
measurements  by  Keq. 


For  this  problem,  the  choice  of  Aj  was  made 
by  applying  eigenstructure  synthesis  to  a  central 
design  flight  condition  and  using  the  resulting  closed 
loop  A  matrix  as  A(j.  The  design  point  used  for  this 
paper  is  trimmed  flight  at  20000  ft  and  a  Mach  number 
of  0.5.  The  resulting  angle  of  attack  is  six  degrees  and 
the  corresponding  A  matrix  is  shown  in  (5.7)  and  the  B 
matrix  in  (4.9). 


OUTER  LOOP  DESIGN 

The  outer  loop  compensator  is  designed  to 
achieve  closed  loop  stability  and  acceptable  performance 
in  the  presence  of  uncertainty.  These  robust  control 
objectives  are  met  by  using  an  implicit  model 
following  formulation  of  p-synthesis.  Figure  3  shows 
the  design  model  used. 
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Figure  3  Design  Model 


The  exogenous  input(w)  is  the  command  from  the 
prefilter,  and  the  feedback  is  the  measured  error  in  p  and 
stability  axis  roll  rate.  The  outputs(e)  consist  of  a 
dynamically  weighted  error  between  ideal  and  actual 
response,  statically  weighted  actuator,  and  dynamically 
weighted  input  command.  DK-iteration"^  is  used  to 
minimize  the  singular  values  of  the  transfer  function. 
The  closed  inner  loop  plant  equalization  ensures  this 
same  compensator  will  achieve  design  goals  for  all 
flight  conditions,  provided  the  inversion  error  is 
sufficiently  small.  The  ideal  model  of  Figure  3  is 


Ideal= 


r  ^  1 

r  3  “I 

^com 

—  0 
s+3 

0  ^ 

_^coin  _ 

s+3_ 

(6.1) 


sin  a 

cos  O'*] 

(6.4) 

-cos  a 

sin  a  J 

[’  ° 

'1 

(6.5) 

[0  cos  a 

sinocj 

The  resulting  controller  is  stable  and  23rd 
order.  This  is  obviously  not  practical  for  application  so 
it  is  reduced  to  10th  order  using  balanced  truncation^ 
Any  reduction  beyond  10th  order  produces  undesirable 
performance  and  robustness. 

COMMAND  SHAPING  AND 
IMPLEMENTATION 

The  desired  handling  qualities  are  achieved 
through  the  use  of  prefilters.  One  can  directly  design 
handling  qualities  into  the  outer  loop  p 
controller^’^"^’^^,  but  having  changing  handling 
qualities  with  different  flight  conditions  would  require 
multiple  scheduled  outer  loop  controllers  which  is 
something  that  we  were  trying  to  avoid  in  the  first 
place.  With  the  p-synthesis  controller,  we  can  assume 
that  the  aircraft  will  respond  with  the  dynamics  shown 
in  (6.1).  This  assumption  allows  us  to  produce  desired 
responses  through  the  use  of  the  prefilter  shown  in 
(7.1). 


Given  a  good  p-synthesis  design,  the  sideslip  and 
stability  axis  roll  rate  responses  will  be  first  order  and 
decoupled  due  to  the  diagonal  nature  of  (6.1).  By 
forcing  the  complementary  sensitivity  function  to  take 
the  frequency  response  shape  of  this  ideal  model,  the 
prefilter  can  provide  the  necessary  handling  qualities  by 
canceling  out  those  dynamics  and  shaping  the  response. 
The  performance  weight  is  chosen  to  provide  close 
following  to  the  ideal  frequency  response  below  ten 
radians/second  and  is: 


Wp 


~0.1(s+50Q) 

s+0.05 

0 


0 

0.1(s+500) 
s+0.05  _ 


(6.2) 


Uncertainty  is  modeled  as  a  complex  multiplicative 
perturbation  at  the  plant  input  with  weighting  function: 


~10(s+100) 

s+10000 

0 


0 

10(s+100) 
s+10000  _ 


(6.3) 


The  actuator  weighting  is  0.0001  time  a  two  by  two 
identity  matrix.  The  matrices  T}  and  T2  (Fig.  3)  are 
transformations.  T2  transforms  p,  p,  and  r  to  (i  and 
stability  axis  roll  rate  while  Ti  transforms  the  outputs 
of  the  p  controller  to  p  and  f  commands.  They  are 
both  dependent  on  angle  of  attack  and  are  shown  in  (6.3) 
and  (6.4). 


flj^(s+3) 

) 


Ideally,  the  zeros  in  the  prefilter  will  cancel  out  the 
aircraft  dynamics  leaving  the  desired  handling  qualities 
response.  Even  if  the  cancellation  is  off,  the  response 
should  be  well  behaved.  Some  care  must  be  taken  in 
the  implementation  of  the  stability  axis  roll  rate 
filtering  since  it  is  a  proper  transfer  function.  As 
mentioned  previously,  it  is  desired  to  change  the 
handling  qualities  of  the  aircraft  with  angle  of  attack. 
This  is  accomplished  by  scheduling  the  parameters  Tj, 
(0,  and  ^  in  (7.1)  with  angle  of  attack.  The  functions 
for  the  three  parameters  are  shown  in  (7.2),  (7.3),  and 
(7.4). 


^com  IJ 
^comj  ! 


3(s^  +2^m+o)^ 


Tj.  =0.0001a^  +0.01 17a+0.3  (7.2) 

6)=-0.04a+3.0  with  (y>0.5  (7.3) 

f=0.005a+0.7  (7.4) 

The  implementation  of  this  control  law  is 

shown  in  Figure  4. 
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Pcmdl 


Ti(a)K^(s) 


,  _  5, 

hO^T(M,h)  H-sfWl 

r  iSvTv' 

Keq(M,h)  ' 


T9(a) 


G(M,h,s)h 


Controller 

Figure  4  Controller  Implementation 


The  pilot  commands  are  fed  through  the  prefilter  and  the 
difference  is  taken  with  feedback  measurements  of 
sideslip  and  stability  axis  roll  rate.  The  error  is  sent 
through  the  |LL-synthesis  controller  to  produce  the  outer 
loop  commands.  The  inner  loop  equalization  takes  the 
sideslip  and  body  axis  rates  and  produces  the  inner  loop 
commands  which  are  then  combined  with  the  outer  loop 
commands.  The  control  selector  transforms  the 
generalized  commands  into  the  actuator  commands.  The 
equalization  gains,  K^q,  and  the  control  selector,  T,  are 
functions  of  Mach,  altitude,  and  a.  The  outer  loop 
controller,  K^,  is  a  single  fixed  dynamic  compensator. 


ROBUSTNESS  ANALYSIS 

In  this  section  a  model  (Figure  5)  reflecting 
structured  and  unstructured  uncertainty  stability 
robustness  requirements  is  evaluated.  The  M-D  block 
diagram  has  inputs  and  outputs  for  a  structured 
uncertainty  block  (wj  -  zj),  an  unstructured  uncertainty 
block  at  the  input  to  the  actuators  (W2  ’  zj),  and  an 
unstructured  uncertainty  block  at  the  output  of  the 
sensors  (W3  -  Z3).  All  perturbations  are  considered  to 
act  simultaneously.  Notice  that  the  pre filter  is  left  out 
of  the  formulation  since  it  is  always  stable  and  will 
produce  stable  inputs  to  the  system  in  Figure  5  for 
stable  pilot  commands. 


The  structured  uncertainty  is  based  on  the 
accuracy  of  the  aerodynamic  database  and  shows  up  in 
the  plant  A  and  B  matrices.  The  level  of  uncertainty  for 
each  derivative  is  shown  in  Table  II  and  is  captured  as  a 
percentage  of  its  nominal  value^. 


Table  II  Structured  Uncertainty  Levels 


Stability  Derivatives 

Control  Derivatives 

Yb=15% 

Y5r=15% 

L3=10% 

^5vtv— 15% 

Ld=30% 

L5de=15% 

Lt=20% 

L5df=10% 

N3=30% 

L5i=40% 

No=50% 

L§  V  t  v~^^  % 

Nr=15% 

N5de=15% 

N5df=20% 

N5r=15% 

N§vtv“15% 

These  uncertainties  are  incorporated  into  the 
analysis  through  the  weights^^ 


BaH 


'1  000000000010010  o' 
01010101010010010 

0  0  1  0  0  1 


0 

0 

0 

0 

0 

0.20L^ 


Ca  = 


0.I5L. 


1010 

1  0 

0101 

0  1 

■0.15Y^ 

O.IOL^ 

0.30N^ 

0 

0 

0 

0 

0.30Lp 

0 

0.50Np 

0 

0 

0 

0 

- 

®10x3 

®7x4 

(8.1) 


0.15N^ 


DE 


0 

O.lOLc 


0.20Nc 


'’DF 

^DF 


0.15N^ 


(8.2) 


0 

0 

O.ISYx 

0 

0 

0 

0.40Lx 

0 

0 

0 

0.I5Njt 

^R 

0 

0 

0 

0 

0.15Y;f 

Oytv 

0 

0 

0 

0.40L. 

^YTV 

0 

0 

0 

0.15No 

^YTV_ 

(8.3) 


Unstructured  uncertainty  at  the  actuator  inputs  is 
assumed  to  be  complex  and  is  weighted  by  the  gain 
matrix 


w 


act 


“0.15 

0 

0 

0 


0 

0.15 

0 

0 


0  0  “ 

0  0 

0.15  0 

0  0.15 


(8.4) 


which  bounds  the  uncertainty  to  15%  over  all 
frequencies.  The  sensor  uncertainties  are  also  considered 
as  complex,  unstructured,  and  weighted  by  the  gain 
matrix 
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(8.5) 


0  0  " 

O.OI  0 
0  0.01_ 

This  weighting  assumes  a  bounded  error  of  10%  for  p 
measurements  and  1%  for  the  rate  measurements. 

Stability  robustness  is  tested  for  six  flight 
conditions  at  high  angle  of  attack.  Table  III  shows  the 
Mach,  altitude,  angle  of  attack,  and  dynamic  pressure  for 
each  of  the  flight  conditions.  The  bounds  of  |Li  to  the 
previously  mentioned  complex  and  real  uncertainties  are 
shown  in  Figure  6. 


Figure  6  |Li  Bounds  for  Stability  Robustness 

The  [i  bounds  seem  to  peak  between  6  and  8  rad/s  with 
the  upper  bound  having  a  magnitude  of  1.15.  Since 
this  is  greater  than  one,  we  are  not  guaranteed  robust 
stability.  Considering  the  fact  that  high  angle  of  attack 
is  inherently  less  stable,  and  that  we  are  considering 
many  uncertainties  simultaneously,  we  can  say  that  the 
system  is  not  particularly  sensitive  since  the  bound  is 
only  slightly  greater  than  one. 


Table  III  Analysis  Flight  Conditions 


Test 

Point 

Mach 

Altitude 

(ft) 

a  (deg) 

q(psf) 

1 

0.2 

10,000 

28.0 

40.8 

2 

0.25 

10,000 

16.8 

63.7 

3 

0.18 

10,000 

35.8 

33.0 

4 

0.25 

20,000 

26.8 

42.6 

5 

0.3 

20,000 

17.4 

61.3 

6 

0.35 

30,000 

20.3 

53.9 

Figure  7  LOES  Dutch  Roll  Frequency  and  Damping 


Analysis  in  the  time  domain  is  accomplished 
through  SIMULINK  simulation  in  MATLAB,  and  the 
only  nonlinearities  considered  here  are  actuator  rate  and 
position  limits.  Since  one  of  the  goals  of  this  design  is 
full  envelope  control,  two  different  flight  conditions  are 
used  as  test  points.  The  first  is  trimmed  flight  at 
M=0.7  at  20000  ft  and  the  second  is  trimmed  flight  at 
M=0.18  at  10000  ft  with  three  degrees  and  thirty  five 
degrees  angle  of  attack  respectively.  The  first  set  of 
responses  are  to  a  1  deg  sideslip  command  and  are 
shown  in  Figures  9,  10,  11,  and  12. 


0.1 

0 

0 


FLYING  QUALITIES  EVALUATION  AND 
SIMULATION  RESULTS 

The  flying  qualities  are  evaluated  by  fitting  a 
LOES  to  the  full  order  closed  loop  transfer  function  as 
discussed  in  Section  III  for  the  6  flight  conditions  given 
in  Table  III.  Figures  7  and  8  show  the  evaluation 
results  graphically,  indicating  that  Level  1  flying 
qualities  are  achieved  at  all  flight  conditions  tested. 
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DegrcMandDegrees/scc  Degr^s  and  Degrees/sec 


A0A=3  degrees,  M=().7,  AIUtudc=200{X)ft 


Figure  9  Sideslip  and  Stability  Axis  Roll  Rate 
Response  to  Ideg  Sideslip  Command 


Figure  11  Actuator  Response  to  Ideg  Sideslip 
Command  at  Low  AOA 


AOA=35  degrees,  M=().l  8,  Allitudc=l()00()ft 


Figure  10  Sideslip  and  Stability  Axis  Roll  Rate 
Response  to  Ideg  Sideslip  Command 


Figure  12  Actuator  Response  to  Ideg  Sideslip 
Command  at  High  AOA 

In  the  responses,  solid  lines  represent  sideslip  and  the 
dash-dot  lines  represent  stability  axis  roll  rate.  The 
ideal  responses  are  shown  as  dashed  lines.  Both 
responses  show  that  the  actual  output  follows  the 
desired  response  quite  closely.  Also,  both  cases  have 
good  decoupling  between  roll  and  sideslip  although,  as 
one  would  expect,  there  is  more  coupling  at  higher 
angles  of  attack.  There  is  also  a  little  bit  of  time  delay 
in  both  cases  which  can  be  attributed  to  the  fourth  order 
dynamics  of  the  actuators.  One  can  also  see  how  the 
prefilter  commands  slower,  higher  damped  responses  at 
higher  angles  of  attack.  Figures  11  and  12  show  how 
the  controller  uses  the  various  effectors  to  accomplish 
the  responses.  Although  the  command  is  the  same  for 
both  cases,  the  control  usage  is  different.  The  higher 
angle  of  attack  condition  requires  more  body  axis  roll 
control  even  though  the  command  is  less  aggressive 
than  the  low  angle  of  attack  condition. 
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Figures  13,  14,  15,  and  16  also  show 
responses  for  the  two  flight  conditions,  but  these  are  to 
a  10  deg/sec  stability  axis  roll  rate  command. 


AOA=3  degrees,  M=0.7,  Altitude=20000ft 


Figure  13  Sideslip  and  Stability  Axis  Roll  Rate 
Response  to  lOdeg/sec  Roll  Rate 
Command 


AOA=35  degrees,  M=().18,  AUiiude=l(X){K)n 


Figure  14  Sideslip  and  Stability  Axis  Roll  Rate 
Response  to  lOdeg/sec  Roll  Rate 
Command 


Figure  15  Actuator  Response  to  lOdeg/sec  Roll 
Rate  Commandat  Low  AOA 


Figure  16  Actuator  Response  to  lOdeg/sec  Roll 
Rate  Commandat  High  AOA 

The  results  for  this  input  are  very  similar  to  the  sideslip 
command  results.  The  response  follows  the  desired 
response  very  well  at  low  angles  of  attack  and 
satisfactorily  at  high  angles  of  attack.  Again,  the 
sideslip  and  roll  are  very  well  decoupled.  The  biggest 
difference  in  the  high  and  low  angle  of  attack  cases  can 
be  seen  in  Figures  15  and  16.  Much  larger  deflections 
are  necessary  for  the  high  angle  of  attack  condition 
which  is  to  be  expected.  Also,  it  can  be  expected  that 
larger  step  inputs  for  the  high  angle  of  attack  case  will 
produce  less  desirable  responses,  given  the  rate  and 
position  limit  information,.  In  all  cases,  it  can  be  seen 
that  the  differential  stabilator  deflection  is  roughly  0.25 
times  the  deflection  of  the  flaperon  and  that  the 
deflection  of  the  yaw  nozzle  is  roughly  0.5  times  the 
rudder  deflection. 
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CONCLUSIONS 

A  robust  flight  control  design  for  the  lateral- 
directional  dynamics  of  the  VISTA/MATV  F-16  has 
been  presented  which  combines  linear  dynamic 
inversion  and  |i-synthesis.  An  inner/outer  loop 
structure  and  pseudo-controls  are  used  to  achieve 
performance  and  robustness  goals  for  a  high  angle  of 
attack  flight  envelope  without  gain  scheduling.  The 
approach  uses  a  command  prefilter,  control  selector,  an 
inner  equalization  loop,  and  an  outer  robust  performance 
loop  to  achieve  excellent  performance  and  robustness. 
Simulation  results  showed  that  good  performance  can  be 
achieved  at  both  low  and  high  angles  of  attack. 
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